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Nitrogen bearing species are common tracers of the physical conditions in a wide variety of objects, and most remarkably in dark 
clouds. The reservoir of gaseous nitrogen is expected to be atomic or molecular, but none of the two species are observable in the 
dark gas. Their abundances therefore derive indirectly from those of N-bearing species through chemical modelling. The recent years 
have accumulated data which stress our incomplete understanding of the nitrogen chemistry in dark cloud conditions. To tackle this 
problem of the nitrogen chemistry in cold gas, we have revised the formation of nitrogen hydrides, which is initiated by the key 
reaction + H, — > NH^ + H. We propose a new rate for this reaction which depends on the ortho-to-para ratio of H2. This new 
rate allows to reproduce the abundance ratios of the three nitrogen hydrides, NH, NHj, and NH3, observed towards IRAS 16293-2422, 
provided that the channel leading to NH from the dissociative recombination of NjH^ is not closed at low temperature. The ortho- 
to-para ratio of H2 is constrained to O/P = 10"^ by the abundance ratio NH:NH2, which provides a new method to measure 0/P. 
This work stresses the need for reaction rates at the low temperatures of dark clouds, and for branching ratios of critical dissociative 
recombination reactions. 
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1. Introduction 

Chemistry in an astrophysical context is not only a matter in 
itself but also provides invaluable tools to determine physical 
conditions such as volume density and kinetic temperature (?). 
Astrochemistry relies on chemical networks which depend on 
the type of environment (e.g. diffuse vs dense gas). The kinetic 
rates of the reactions are, in the best cases, based on thermody- 
namical and quantum mechanical calculations and experiments, 
which make the setup of such chemical networks an extremely 
demanding process. Various such networks of reactions are pub- 
licly available (KIDA, UMIST, OSUQ, Flower & Pineau des 
Foretfl), with various degrees of complexity and/or complete- 
ness regarding specific aspects of the chemistry (e.g. deuteration, 
cations...). 

The nitrogen element is among the 5th or 6th most abun- 
dant in the Solar neighbourhood, after H, He, C, O, and probably 
Ne (??). In the cold neutral medium, nitrogen is expected to be 
predominantly atomic or molecular, but direct observations of 
N or Nj are not possible. The amount of gaseous nitrogen thus 
relies on the abundances of N-bearing molecules via chemical 
models. Nitrogen bearing molecules are observed in a wide va- 
riety of physical conditions. Molecules such as CN, HCN, and 
HNC, with large permanent dipole moments, are detected to- 
wards diffuse clouds (?), dense cores (?), protoplanetary disks 
(?) and high-z galaxies (?). N-bearing species, such as NjH^ 
(and its deuterated isotopologues), are also efficient tracers of 
the dense and cold gas where CO has already frozen-out (???). 
The cyanide radical CN is also a precious molecule which serves 
as a tracer of the magnetic fields through Zeeman splitting (?). 
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Understanding the chemistry of nitrogen is thus crucial in many 
astrophysical areas. 

Ammonia and NjH^ are daughter molecules of N,, which 
forms from atomic N through neutral-neutral reactions medi- 
ated by CN and NO (?). The formation of CN and the related 
HCN and HNC molecules, and in particular their abundance ra- 
tio CN:HCN, are however not fully understood (?). The forma- 
tion of ammonia in dark and dense gas is thought to take place in 
the gas phase. Once Nj exists in the gas phase, it reacts with He^ 
ions formed by cosmic-ray ionization, to form N^ which by sub- 
sequent hydrogen abstractions lead to NH^ and finally NH3 by 
dissociative recombination. ? (LB91 in the remainder of the pa- 
per) was the first to investigate the influence of the Hj ortho-to- 
para ratio (noted O/P) on the formation rate of NH3. LB91 con- 
sidered the conversion between p-Hj and o-Hj through proton 
exchange in the gas phase. He concluded that when O/P > 10"^, 
the formation of ammonia no longer depends on the value of 
O/P. However, a better test of the nitrogen chemistry is pro- 
vided by the simultaneous observations of all three nitrogen hy- 
drides, NH, NHj, and NH3 (e.g., ?). The Herschel/HIFI instru- 
ment has opened the THz window which contains the funda- 
mental rotational transitions of light molecules with large elec- 
tric dipole moments such as hydrides. ?, Paper I in what fol- 
lows have detected NH, NHj, and NH3 in the cold envelope of 
the Class protostar IRAS 16293-2422. The hyperfine structures 
of the = 1 - transitions of NH and NHj are seen for the 
first time, in absorption, and allow a precise determination of 
their excitation temperature (Tex) and line centre opacity. Four 
lines of ammonia are also detected in absorption, plus the funda- 
mental rotational transition at 572 GHz which evidences a self- 
absorbed line profile, likely tracing both warm and cold gas. For 
all three molecules, excitation temperatures were found in the 
aiah|;e 8-10 K (Paper I). The derived column density ratios are 
NH : NH2 : NH3 =5:1: 300. Whilst abundances are generally 
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Fig. 1. Comparison of the rate coefficient as given by Eq.[3]to ex- 
perimental data for the reaction of with n-H, and p-H, with 
different o-Hj admixtures: 13%, 3%, 0.8%, andO% (blue'hne). 
The small symbols are taken from ? while the large triangles are 
CRESU results from ?. The solid lines correspond to Eq. [3]of 
the present paper. 

delicate to derive because the column density of is uncertain, 
the column density ratios above provide a stringent test for the 
first steps of the nitrogen chemistry and for the ammonia syn- 
thesis in particular. Indeed, the current nitrogen networks avail- 
able all produce more NHj than NH under dark cloud conditions 
(Paper I), at odds with the observations. 

In this Letter, we investigate the NHiNHj problem in dark 
gas, by revisiting the kinetic rate for the reaction 

N+-HH2 — >NH+-i-H (1) 

We derive separate rates for reaction with p- and o-Hj for which 
the nuclear spin 7 = and 1 respectively. In Section 2 we de- 
scribe the new rates which are compared to experimental mea- 
surements and derive the abundances of nitrogen hydrides in typ- 
ical dark cloud conditions. The results are discussed in Section 3 
and we propose conclusion remarks in Section 4. 

2. Chemical modelling 

2.1. Rate of reaction + 

Reaction (1) has a small activation energy in the range 130- 
380 K (?), with a consensus now tending towards a value 
below 200 K (e.g. ?). This reaction has been studied ex- 
perimentally by ? who considered pure para-Ha and a 3:1 
mixture of o- and p-Hj, refered to as normal H2 (n-Hj). 
These authors fitted the two rates with normal and P-H2 
as k„ = 4.16x10""' exp[-41.9/r] cm^s"' and = 
8.35x10""' exp[-168.5/r] cm^s"'. The rate with p-Hz 
was used by LB91, who assumed that the reaction with o-Hj 
proceeds with no endothermicity, on the basis that the 170.51 K 
internal energy of the J - 1 level of o-Hj is used to overcome 
the endothermicity. Accordingly, the rate of reaction ([TJ was 
written as /tLB9 1 = 8.35 xlO""'[jt;-H(l -x) exp(-168.5/r)], where 
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Fig. 2. Comparison of the rate coefficient for reaction ((TJ as 
given by Eq. [3] to the rates from LB91 (dashed) and the OSU 
database (dot-dashed). The rates are computed for temperatures 
in the range 8-150 K and O/P = 10 ^ 10 ^ 0.01, 0.1, and 3 
(black lines). The symbols show the fits from ? for normal H2 
(red circles, 0/P=3) and P-H2 (blue triangles). 

X - 0/P/(l + O/P) is the fraction of o-Hj. Actually, the expres- 
sions for k„ and kp allow to derive a simple expression for ko, 
the reaction rate of with o-Hj, namely k„ - 3/4ko + l/4-kp 
(?). For 0/P=3:l and at temperatures less than w 40 K, the rate 
is essentially k,,, and k„ is thus a good approximation of k„ at 
the low temperatures (T < 15 K) of dark clouds. Alternatively, a 
single-exponential fit to kg leads to (see Fig. lA. Il l 

A:o = 4.2x10"'° (r/300)"*'-'^ exp[-44.5/r] cm^ s"' (2) 

The rate for reaction ([T]i with an H2 admixture of arbitrary O/P 
is then obtained as 

ki - xkg + (I - x)kp (3) 

In Figure[T] the rate ki from Eq.[3]is compared to the experimen- 
tal data of ? for the reaction with n- and p- Hj down to 8 K and 
45 K respectively. ? performed experiments with n-Hj and with 
p-H, containing admixtures of o-Hj (13%, 8%, and 0.8%), at 
temperatures down to 14 K. The agreement between our Eq. [3] 
and the three sets of p-Hj data of ? is excellent at temperatures 
below 20 K, and to within a factor of 2 up to 100 K. This suggests 
that Eq. [3] should be accurate to within a factor of 2-3 down to 
10 K. On the other hand, this rate should not be employed above 
~ 150 K. 

We compare, in Fig. |2] the rate given by Eq. [3] to 
fcLB9i and to the rate given by the OSU database A;osu = 
8.35 xlO""' exp[-85/r] cm"* s"'. The rate is computed for sev- 
eral values of O/P. As expected, it converges towards the n-Ha 
and P-H2 rates for large (0/P=3) and small (0/P= 10"'') val- 
ues respectively. It is evident that LB91 overestimated the rate 
of reaction (1) by several orders of magnitude at temperatures 
smaller than 20 K. In warmer gas (T > 20 K), k[ and fcLB9i agree 
to within a factor of 10 or less, whilst the OSU rate approaches 
the 0/P=3 curve. The OSU rate thus amounts to assume O/P 
ratios larger than 0.01, whereas LB91 provides acceptable rates 
only at temperatures larger than 20 K. 
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Table 1. Initial gas phase fractional abundances (n(X)/n}i) from 
? with «H = n(H) + 2n(H2). For comparison, the fractional abun- 
dances from the low-metal model of ? are also given. Numbers 
in parentheses are powers of 10. 



Species 


Flower 2003 


Wakelam 2008 


H2 


0.50 


0.50 


He 


0.10 


0.14 


N 


6.39(-5) 


2.14(-5) 





1.24(-4) 


1.76(-4) 


c+ 


8.27(-5) 


7.30(-5) 


s+ 


1.47(-5) 


8.00(-8) 


Si+ 


8.00(-9) 


8.00(-9) 


Fe+ 


3.00(-9) 


3.00(-9) 


Mg+ 


7.00(-9) 


7.00(-9) 



2.2. Abundance of nitrogen hydrides 

In this work, we wish to reproduce the ratios NHiNHjiNHj ob- 
served in the cold envelope of IRAS 16293-2422 (Paper I). The 
fundamental hyperfine transitions of NH, NHj and several tran- 
sitions of NH3 have been detected in absorption. Absorption is 
interpreted as resulting from the low temperature of the gas in the 
envelope seen against the warmer continuum emitted by the dust 
closer to the protostar. The gas density is low ( 10**- 10^ cm""*) 
as compared to the critical density of the detected transitions 
(> 10^ cm"^) which ensures that collisions do not govern the 
(de)excitation processes. The = 1 - transitions are there- 
fore thermalized with the dust emission temperature which, in 
the THz domain, correspond to excitation temperatures close to 
a kinetic temperature of 10 K. The lines are Gaussian and do 
not show signatures of strong dynamical effects such as infall. 
Dynamical timescales are then expected to be large with respect 
to the free-fall time. The observed lines therefore trace a cold 
gas, moderately dense, free of dissociating photons, and where 
the ionization is driven by cosmic rays. 

The main point raised by Paper I is that the observed ra- 
tios NH:NH2:NH3=5:1:300 could not be reproduced by chemi- 
cal networks updated regarding the rates of the dissociative re- 
combination (DR) reactions leading to NH, NH,, and ammonia. 
Paper I considered three models, with varying branching ratios 
for some DR reactions. Whilst the NH2:NH3 abundance ratio 
could be reproduced in all three cases, no model was able to 
product! [NH] > [NHj]. These models used the OSU rate for 
reaction ([T]i which was shown in Sec. 2. 1 to depart from the mea- 
sured rate by several orders of magnitude in cold gas and O/P 
smaller than 0.01. The next Section explores the consequence of 
the new rate given by Eq. [3] on the abundances of the nitrogen 
hydrides, in typical dark cloud conditions. 

We have performed chemical calculations in a gas at T = 
10 K, with density hh = 10^ cm""*. The gas is screened by 
10 mag of visual extinction, such that the ionization is primarily 
due to cosmic rays at an adopted rate f = 1.3 xlO"'^ s"'. Gas- 
phase abundances are computed as a function of time by solving 
the chemical network until steady-state is reached. Freeze-out of 
gas-phase species onto dust grains are ignored in this work. In 
what follows, the quoted abundances correspond to the steady- 
state. The initial fractional elemental abundances (nx/ne) are the 
gas phase abundances taken from ?, hereafter FPdF03 who con- 
sidered the depletion of metals in grain mantles, grain cores, and 

^ Hereafter, the fractional abundance of species X is noted [X] and 
equals n(X)/«H- 



PAHs (see Table[TJ. These abundances differ from those adopted 
in Paper I (?). The rates for the dissociative recombination reac- 
tions are those of the model 1 of Paper I. From the above, the rate 
for the key reaction O depends on the O/P which is not known, 
but which might strongly differ from the thermodynamical equi- 
librium value (???) which is ^ 10"' at 10 K. The steady-state 
abundances of NH, NHj, and NH3 have thus been computed for 
various values of O/P in the range 10"' to 3. 

In this work, the DR of NH| has three output channels, 
NH3 + H, NH2 + Hj, and NHj -H H H- H (?). However, another 
channel may be NH -H H H- Hj (?), but to our knowledge, no 
branching ratio is available in the litterature. Values up to 10% 
may be considered in a future work. 

We first consider the case where the DR of NjH^ pro- 
duces only N, + H (?). The resulting abundances are the dashed 
lines in Fig. [3] It appears that the O/P controls the abundances 
of NH, NHj, and NH3 but the ratios NH:NH2 and NH3:NH2 
are insensitive to O/P. In a second series of calculations, the 
NH + N channel is given a 10% branching ratio for the DR of 
NjH^. In this case, the abundance of NH is nearly independent 
of O/P, while the abundances of NHj and NH3 remain unaf- 
fected (full lines in Fig. |3] As a consequence, it is found that 
for O/P < 0.03, NH : NHj > 1. In Paper I though, opening 
this channel did not solve the NHiNHj problem, which contrasts 
with the above result. This may be understood as follows. When 
the NH -I- H channel of the dissociative recombination of NjH^ 
is opened with a 10% braching ratio, the corresponding rate is 
w 5.5 xlO"^ cm-^ s"', several orders of magnitude larger than ki 
(see Fig. |2]i. Hence this channel, which is insensitive to O/P, 
dominates the formation of NH over NHj + e" . In contrast, the 
other hydrides NHj and NH3 are daughter molecules of NH^ 
which is formed from N^ + Hj whose rate does explicitely de- 
pend on O/P in the present work. Hence the abundances of NHj 
and NH3 do also depend on O/P. On the other hand, at the tem- 
perature of 10 K and O/P < 10""^, the revised rate ki (Eq.O is 
at least an order of magnitude smaller than the OSU rate used in 
Paper I. For larger O/P, the rate is similar or even larger up to 
a factor of 10. These two effects, namely the drop of ki for low 
O/P, and the O/P-independant formation of NH, makes it possi- 
ble to produce a NHiNHj ratio with values greater than unity. In 
the following we keep the NH + N channel opened to 10% and 
explore some consequences of this result. 



3. Discussion 

3.1. Tlie H2 O/P ratio in dark clouds 

Computing the abundance ratios NHiNHj and NH3:NH2 for dif- 
ferent O/P ratios leads to the left panel of Fig. |4] The initial frac- 
tional abundances are kept fixed, while the rate of reaction ([T]i is 
calculated for values of O/P from 10"' (close to the Boltzmann 
value at 10 K) to 3. The calculated ratios are compared with 
the observational constraints. As expected, the NH3:NH2 ratio 
is insensitive to O/P variations. It is equal to 10, not consistent 
with the observed value of 300. On the contrary, the NH:NH2 
ratio shows two regimes: for O/P < 10"^, NH : NHj ~ 50, 
and for O/P > 0.1, NH : NHj ~ 0.07. The observed ratio of 
5 is intermediate, and indeed well constrains the O/P to 10"^. 
From Fig. |4] we conclude that there is a range of O/P for which 
[NH] > [NHj] and [NH] < [NH3] simultaneously. In our case, 
this range is O/P = 5 xlO"'* - 0.02. The O/P ratio appears as a 
control parameter for the NH:NH2 ratio, and allows for the first 
time to produce more NH than NHj- Conversely, the measure of 
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Fig. 3. Steady-state abundances (with respect to H nuclei) of 
nitrogen hydrides as a function of O/P, in a 10 K gas with 
nH= lO'* cm"^, and ^ - 1.3x10 s"'. The rate of reaction 
(1) is given by Eq. [3] Two branching ratios for the channel 
NjH^ -He" — > NH -H N are considered: 0% (dashed lines) and 
10% (full lines). The observed abundances from Paper I are also 
indicated (filled rectangles). 

NHiNHj tightly constrains the O/P ratio. We note that the frac- 
tional abundance of NH3 of 1.9 xlO measured by ? towards 
the starless core L 1544, corresponds to an O/P - 10"^ - 10"^. 
This ratio is a factor of 10 larger than the corresponding ratio 
from the model of LB91. It is thus found that reaction ([T]i reg- 
ulates the formation of NHj and NH3 whilst the formation of 
NH is controlled by the dissociative recombination of NjH^ pro- 
vided the branching ratio towards NH is 10%. The value of this 
branching ratio is uncertain but is likely non-zero (?). 

3.2. Dependence on the initial abundances 

The NHjiNHj ratio is w 10 in the above models, a factor 30 
below the observed values towards IRAS 16293-2422. However, 
models computed in Paper I with the initial abundances of ? lead 
to a ratio close to the observed value, hence suggesting that the 
initial elemental abundances influence this ratio. The gas phase 
elemental abundances of C, N, and O in dark clouds are poorly 
known because the amount of these elements incorporated into 
the dust (core or mantles) is loosely constrained. The elemental 
abundance of oxygen in the gas phase is not known accurately, 
and variations by an order of magnitude are fully conceivable 
(?), whilst that of carbon is better known. Accordingly, we have 
considered variations of the initial elemental abundances of oxy- 
gen, to vary the ratio C:0, encompassing the value of 0.41 from 
7 

The resulting abundance ratios, computed for an O/P ratio 
of 10"^, are shown in Fig. |4] (right panel). As C:0 decreases 
below the 0.66 value of FPdF03, the ratio NHjiNHj increases. 
In the process, NHiNHj remains constant. When C:0 is now 
increased above 0.66, both ratios decrease, by at most a factor 
2 to 5. It is thus apparent that, when [O] is increased, the C:0 
ratio controls the NH2:NH3 ratio. On chemical grounds, the C:0 
ratio is expected to influence the abundance of NHj, for which 
the main destruction routes involve oxygen to form NO, NH, and 
HNO. Similarly, NH is mostly removed by reaction with oxygen 
to form principally NO. However, when oxygen is significantly 
depleted from the gas phase (C:0 > 0.5), another destruction 



route of NH, involving sulfur, becomes important. As a result, 
the NHiNHj ratio is mostly insensitive to C:0, until C:0> 0.5 
when it starts to decrease by small factors. 

The situation is different for NH3 which is destroyed by 
charge transfer reactions with H^ and S^ and by proton ex- 
change reactions with H3 and HCO^, forming notably NH|, 
which dissociates back into NH3, and NH3 . A fraction of NH3 
will lead to NH, and NH which are the true destruction chan- 
nels of ammonia. The abundance of NH3 is thus only marginally 
affected by the change in C:0, unless the ionization fraction 
is modified, which is the case for low C:0 ratios. An increase 
of the oxygen abundance (i.e. a decrease of C:0) is accompa- 
nied by a deacrease of [e^], or equivalently [S^] - the dom- 
inant ion - which makes the abundance of NH3 to increase. 
Consequently, as C:0 decreases below 0.66, NHiNHj keeps con- 
stant and NH3:NH2 increases by more than one order of magni- 
tude. Whereas when C:0 increases, NHiNHj decreases by less 
than a factor 10, whilst NH3:NH2 is approximately constant. 

4. Conclusions 

Using an updated rate for reaction ^ with an explicit de- 
pendence on the O/P ratio, we have shown that the O/P of 
H2 controls the ratio NHiNHj in dark clouds without affect- 
ing the NH3:NH2 ratio. A value of O/P close to 10"^ leads to 
NH : NHj > 1 and NH3 : NHj < 1 as observed. Interestingly, 
this value of O/P is close to the predictions of ?, their Fig. 1 
under similar conditions. In addition, measuring the NH:NH2 
ratio may be a new method to constrain the O/P ratio of H2 in 
dark clouds. We have also shown that decreasing the C:0 ratio 
by increasing [O] controls the NH2:NH3 ratio. Finally, accept- 
able parameters are found that lead to abundance ratios in agree- 
ment with the observations towards IRAS 16293-2422, namely 
O/P ~ 10"^ and C:0< 0.4. It is to be noted, however, that in 
this range of parameters, the absolute abundances predicted by 
the models are a factor of 10 below those derived in Paper I. 

In this work, the influence of O/P on the ratios NH:NH2:NH3 
has been explored only through the rate of reaction ([TJ. A bet- 
ter study would be to self-consistently compute the O/P ratio 
from a model of conversion of o-Hj into p-Hj, as was first con- 
sidered by LB91, and more recently by e.g. ?. Including proton 
exchange reactions, in the gas phase, between H2 and H^ (?), Hj 
(?), and H3 (?), would reprocess the 3:1 mixture of H2 formed 
on grains to a different O/P. An even more self-consistent ap- 
proach would be that of ? who separate reactions with the ortho- 
and para- forms of all concerned N-bearing molecules. Another 
obvious continuation would be to explore these findings in dif- 
ferent physical conditions (Ay, cosmic -ray ionization rate, etc). 

Branching ratios of dissociative recombination reactions are 
crucial to the formation of nitrogen hydrides. In this work, the 
dissociative recombination of NH4 has three output channels 
NH3 + H, NH2 -I- H2, and NH^ -H H -1- H (?). However, another 
channel may be NH -H H H- Hj (?), but to our knowledge, no 
branching ratio is available in the litterature. Values up to 10% 
may be considered in a future work. Critical too are the values of 
rates at low temperatures, especially for reaction (1) with p-Hj 
for which measurements at temperatures lower than 14 K are not 
available. 
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Fig. 4. Abundance ratios calculated when a 10% branching ratio is considered for the reaction NjH^ + e" — > NH + H. The hashed 
bands show the observed ratios with their Icr uncertainties. Left panel: varying O/P, with the initial abundances of ? ([0] = 1.24(-4) 
and [C]=8.27(-5).)- Right panel: varying [O] while keeping [C] constant, at constant 0/P= 10"^. The low-metal abundances of ? 
are [0]=1.76(-4) and [C]=7.30(-5), or C:O=0.41. 
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Fig.A.l. Rate for reaction ([TJ with o-Hj and fit result (see 
Eq.|2|. Bottom panel: filled squares are the values derived from 
?, and the fitted curve (full line). Top panel: relative error. 

Appendix A: Reaction rate 

For the purpose of implementing the rate of the reaction with ortho-H2 in 
chemical networks. Figure IaTT] shows the result of a single-exponential fit. The 
filled squares are based on the fitted reaction rates from ?. A weighting in the 
form l/(r - 10)^ was adopted, which is accurate to better than 6% for T = 
8 - 150 K. 



